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Abstract. Over the past few years three dimensional Quantum Electrodynamics (QED3) has 
attracted a lot of attention, because it may be an effective theory for the underdoped and non- 
superconducting region of the phase diagram of high T c cuprate compounds. We present results 
from lattice simulations of the non-compact version of the theory in order to address the issue 
of whether chiral symmetry is spontaneously broken when the number of fermion flavours Nf 
is less than a critical value Nf c . Our results provide strong evidence that QED3 is chirally 
symmetric for Nf > 1.5, implying that a pseudogap phase separates the superconducting phase 
from the antiferromagnetic phase. 



1. Introduction 

Interest in QED3 has recently been revived by the suggestion that the model may be an effective 
theory for the underdoped and non-superconducting region of the phase diagram of high-T c 
superconducting cuprate compounds [HE]- In brief, superconductivity in these substances 
is confined to planes defined by Cu02 layers, thus motivating a (2 + l)d description. The 
superconducting order parameter has a d-wave symmetry, implying that there are four nodes in 
the gap function as the Fermi surface (which in (2+l)d is a curve) is circumnavigated. At each 
node the low-energy quasiparticle excitations obey an approximately linear dispersion relation 
with the result that it is possible to rewrite the action for eight distinct low energy species (spin 
up and spin down at each of four nodes) in a relativistically invariant form in terms of Nf = 2 
species of four-component Dirac spinors. For phenomenologically relevant models the action for 
each individual flavor exhibits a spatial anisotropy, a feature ignored in this paper. If QED3 is a 
relevant effective theory for cuprates, then the abstract theoretical problem of the value of Nf c 
assumes concrete phenomenological reference. If Nf c > 2, then the theory is chirally broken at 
zero temperature. On retranslating from the Dirac spinor basis to the original electron degrees 
of freedom, the chiral order parameter is reinterpreted as an order parameter for spin density 
waves, whose wavevector gets shorter and shorter as doping is decreased, until at zero doping 
the Neel antiferromagnetic state is recovered [5J. This picture therefore predicts the existence 
of a phase boundary between superconducting (dSC) and antiferromagnetic (AFM) phases at 
some non-zero doping in the zero temperature limit. If, on the other hand, Nf c < 2, the 
chirally symmetric ground state manifests itself as a tongue of "pseudogap" phase separating 
dSC from AFM, in which normal Fermi liquid properties may be modified as a result of a 



non-perturbative anomalous dimension for the fermion field PQ. In addition, QED3 has found 
interesting applications in the unconventional quantum Hall effect in graphene [3]. 

The study of quantum field theories in which the ground state shows a sensitivity to 
the number of fermion flavours Nf is intrinsically interesting. Apparently, for Nf > Nf c , 
the attactive interaction between a fermion and an antifermion due to photon exchange is 
overwhelmed by the fermion screening of the theory's electric charge. Initial studies of QED3 
based on Schwinger Dyson equations (SDEs) using the photon propagator derived from the 
leading order 1/Nf expansion suggested that for Nf less than Nf c ~ 3.2 chiral symmetry is 
broken [3J. Other studies taking non-trivial vertex corrections into account predicted chiral 
symmetry breaking for arbitrary Nf [5]. Studies which treat the vertex consistently in both 
numerator and denominator of the SDEs have found Nf c < 00, with a value either in agreement 
with the original study [BJ, or slightly higher Nf c ~ 4.3 [7j. An argument based on a 
thermodynamic inequality predicted Nf c < | [8], a result that was later challenged in [9]. 
Progress in the direction of gauge covariant solutions for the propagators of QED3 showed 
that in the Landau gauge a chiral phase transition exists at Nf c ~ 4 [10] . A gauge invariant 
determination of Nf c based on the divergence of the chiral susceptibility gives Nf c ~ 2.16 [llj . 
It has also been shown that the issue of the gauge dependence of Nf c extracted from SDEs 
becomes irrelevant if Landau-Khalatnikov-Fradkin transformations are taken into account [12] . 

Recent lattice simulations showed that chiral symmetry is broken for Nf = 1, whereas Nf = 2 
appeared chirally symmetric with an upper bound of 10~ 4 on the dimensionless condensate [13] . 
The principal obstruction to a definitive answer has been large finite volume effects resulting 
from the presence of a massless photon in the spectrum, which prevent a reliable extrapolation 
to the thermodynamic limit. Recent lattice simulations of the three-dimensional Thirring model, 
which may have the same universal properties as QED3, predicted Nf c = 6.6(1) [H]. In this 
paper we present preliminary results from lattice simulations of QED3 on large lattices in an 
effort to detect chiral symmetry breaking for Nf = 0.5, 2. 

2. Lattice Model and Simulations 

We are considering the four-component formulation of QED3 where the Dirac algebra is 
represented by the 4x4 matrices 70, 71 and 72. This formulation preserves parity and gives 
each spinor a global U(2) symmetry generated by 1, 73,75 and 27375; the full symmetry is then 
U(2Nf). If the fermions acquire dynamical mass the U(2Nf) symmetry is broken spontaneously 
to U(Nf) x U(Nf) and 2A r | Goldstone bosons appear in the particle spectrum. 
The action of the lattice model we study is 

a N 

S = I E &^(x)e f , u (x)+J2T,U^)M(x,x') Xi (x') (1) 

x,fi<u i=l x ,x' 

This describes interactions between ./V flavours of Grassmann-valued staggered fermion fields 
X, X defined on the sites x of a three-dimensional cubic lattice, and real photon fields 8 X/Ji 
defined on the link between nearest neighbour sites x, x + fx. Since G 2 is unbounded 
from above, eq.((T]) defines a non-compact formulation of QED; note however that to ensure 
local gauge invariance the fermion-photon interaction is encoded via the compact connection 
U Xfi = exp(i8 x ^), with U x +n—n = U*^. In the fermion kinetic matrix M the Kawamoto-Smit 
phases = (— l) Xl ~^ ^m- 1 are designed to ensure relativistic covariance in the continuum 

limit, and m is the bare fermion mass. 



If the physical lattice spacing is denoted a, then in the continuum limit ad — > 0, eq.([T]) can 
be shown to be equivalent up to terms of 0(a 2 ) to 

S = Y,$ Wt> + igAJ + mW + jF^F^ (2) 

3=1 

ie. to continuum QED in 2+1 euclidean dimensions, with ip, tp describing Nf flavours of four- 
component Dirac spinor acted on by 4x4 matrices 7-, and Nf = 2N. The continuum photon 
field is related to the lattice field via 6 X [i — QqA^ (x) , with dimensional coupling strength g given 
by g 2 = (a/3) -1 , and the field strength F^ u = d^Ay — d^A u . The continuum limit is thus taken 
when the dimensionless inverse coupling (3 — > 00. 

Our numerical simulations were performed using the standard Hybrid Molecular Dynamics 
algorithm. We checked the effects of lattice discretization on the values of the chiral condensate 
by comparing data extracted from simulations at fixed physical volume (L//3) 3 and fixed physical 
mass (3m [16] . The results show that the lattice discretization effects are small for Nf > 0.5 at 
f3 = 0.90, m = 0.005 on 54 3 lattices, whereas for Nf = 0.5 there is an ~ 8% discrepancy between 
the values of the dimensionless condensate /3 2 (^V) at /3 = 0.90 and f3 = 1.20. 
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Figure 1. (3 2 (^ip) vs (3m for Nf = 1.5. Figure 2. /3 2 (^V) vs P m f° r Nf = 2. 

We fitted f3 2 {$ip) at different values of [3m for Nf = 1.5, N f = 2 and fixed f3 = 0.90 to 

(3 2 (^) =ao + ai -((3m). (3) 

The finite size effects are small especially for Nf = 2 [see Figs. ([T]) and (j2J)]. For Nf = 1.5, 80 3 
we extracted from eq.([3]) ao = —6 x 10~ 7 with a statistical error 8 x 10~ 7 , whereas for Nf = 2, 
80 3 we extracted ao = 1.5 x 10 -6 with a statistical error 10 -6 . These results provide evidence 
that QED3 with Nf = 1.5 and Nf = 2 is chirally symmetric with an accurancy 10~ 6 . 

Next, we fitted the values of the dimensionless condensate at different Nf, m, and lattice sizes 
to a renormalization group inspired equation of state that includes a finite size scaling term [TS] : 

m = A({(3 - (3 C ) + CL^)((3 2 (^)r + B([3 2 {^)) s , (4) 

where p = 5- l/(3 m . The results extracted from this fit are: A = 0.0477(38), B = 0.79(2), C = 
10.7(8), Nf c = 1.52(6), 6 = 1.177(7), p = 0.73(2). The data and the fitting functions are shown 
in Fig. [3l These results are consistent with a second order phase transition. 
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Figure 3. Fits of /3 2 vs. Nf to a finite volume scaling form of the equation of state. 



3. Summary 

The extrapolations of ^('ipip) vs (5m to the chiral limit on lattices with small finite size effects 
show that QED3 with Nf > 1.5 is chirally symmetric with an accuracy of O(10~ 6 ), which is 
in agreement with the theoretical prediction of [8]. This may imply that the Nf = 2 chirally 
symmetric ground state manifests itself in the phase diagram of the cuprate compounds as a 
tongue of "pseudogap" phase separating the superconducting from antiferromagnetic phases, 
in which normal Fermi liquid properties may be modied as a result of a non-perturbative 
anomalous dimension for the fermion field pQ. The preliminary results extracted from fits to 
a finite volume equation of state are consistent with a second order phase transition scenario 
at Nf c ~ 1.5. However, as we mentioned in the previous section lattice discretization artifacts 
are not negligible for Nf < 0.5. This implies that Nf c could be even smaller than 1.5. We are 
performing simulations closer to the continuum limit to clarify this issue. 
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